Magnetic field dependence of hole levels in self-assembled InAs quantum dots 
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Recent magneto-transport experiments of holes in InGaAs quantum dots [D. Reuter, P. Kailuweit, 
A.D. Wieck, U. Zeitler, O. Wibbelhoff, C. Meier, A. Lorke, and J.C. Maan, Phys. Rev. Lett. 94, 
026808 (2005)] are interpreted by employing a multi-band k • p Hamiltonian, which considers the in- 
teraction between heavy hole and light hole subbands explicitely. No need of invoking an incomplete 
energy shell filling is required within this model. The crucial role we ascribe to the heavy hole-light 
hole interaction is further supported by one-band local-spin-density functional calculations, which 
show that Coulomb interactions do not induce any incomplete hole shell filling and therefore cannot 
account for the experimental magnetic field dispersion. 



PACS numbers: 73.21.La,73.22.-f,73.23.Hk 

Knowledge of the electronic structure of semiconduc- 
tor quantum dots (QDs)i is essential to understand and 
predict their physical properties. The energy struc- 
ture of conduction band electrons confined in QDs has 
been thoroughly investigated using capacitance-voltage 
(C-V) 2 and far-infrared absorption spectroscopy. 3 How- 
ever, direct and clear access to the energy structure of 
valence band holes has long been hindered by the high 
density of states arising from their larger effective mass. 4 
Only very recent C-V spectroscopy experiments have 
succeeded in resolving a number of energy levels corre- 
sponding to holes confined in In(Ga)As self-assembled 
QDs. 5,6,7,8 Of particular interest arc the results reported 
by Reuter et al.^£ which studied the energy dispersion 
of hole charging peaks in QDs with a magnetic field B 
applied along the growth direction. It was found that, 
for QDs charged on average with one or two holes, the 
slope of the peaks dispersion was approximately zero. 
When the QDs were charged with three to six holes, the 
peaks shifted towards lower and higher energy in alter- 
nating fashion, the slope of the fifth and sixth peaks be- 
ing about twice as strong as that of the third and fourth 
ones. By using a one-band effective mass Hamiltonian, 
such results were interpreted as a non-sequential shell 
filling, leading to highly spin-polarized ground states* 7 " 
It was speculated that this phenomenon, which has no 
analogy in conduction band electrons, might be due to 
the Coulomb exchange interactions exceeding the rather 
small kinetic energy of holes. 

In this paper, we provide theoretical support for these 
experiments. First, we study the effect of Coulomb corre- 
lations on the electronic configuration of the hole ground 
state in self- assembled In(Ga)As QDs. To this end, 
we use a one-band local-spin-density functional (LSDF) 
model. We find that, even for a rather large QD and pure 
heavy hole effective mass, a sequential filling of the energy 
levels is predicted. It then follows that an explanation for 
the experimental magnetic field-induced dispersion of the 
charging peaks is unattainable within the one-band pic- 



ture. Next, we study the role of heavy hole-light hole 
(HH-LH) coupling, using a four-band Luttinger-Kohn 
Hamiltonian. Our results show that HH-LH interaction, 
together with the spin Zceman splitting, yield a plausi- 
ble interpretation of the experimental data even at an 
independent particle level. 

Structural information of self- assembled QDs is diffi- 
cult to assess, and this imposes a severe restriction for 
a quantitative description of their physics. Nonetheless, 
the sample of QDs studied in RefsQH showed a narrow 
size distribution, and a rough guess of their shape, size 
and composition is possible. Comparing the growth pa- 
rameters and the photoluminescence emission wavelength 
to other experiments and theoretical work, the QDs are 
believed to be lens shaped, the base plane being approx- 
imately circular but slightly elongated along the [0,1,1] 
direction. The QDs height after overgrowth is believed 
to be between 4 nm and 6 nm, and the diameter between 
32 nm and 40 nm. The QDs are not pure InAs but have 
an inhomogeneous InGaAs composition. 9 

In order to see whether Coulomb correlations can in- 
duce an incomplete shell filling in such structures, we in- 
vestigate the electronic configuration and addition energy 
of the V-hole ground state at zero magnetic field (JV=1- 
6) using an effective mass LSDF models The theoretical 
model is fully three-dimensional, which allows us to rep- 
resent the QDs as circularly symmetric lenses with finite 
height, and use bulk material parameters without any 
need to fit them in order to account for the missing verti- 
cal motion.Ai We probe two model QDs, which constitute 
a lower and an upper estimate of the possible experimen- 
tal sizes. The smaller (larger) QD has a radius of R = 16 
(20) nm, and a height of h = 4 (6) nm. In both cases, 
we take an average effective mass value between that of 
InAs and that of GaAs heavy holes, m* = 0.44^ The 
confinement potential is defined by a step-like function 
whose height is given by the InGaAs valence band-offset, 
380 meVji 6 " and the effective dielectric constant we use is 
12.4. 
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Figure 1 represents the addition spectra, AA(N) = 
E(N + 1) - 2 E(N) + E(N - 1), where E(N) is the en- 
ergy of the iV-hole ground state, for the smaller (solid 
line) and larger (dashed line) QDs. The insets show the 
occupation of the Kohn-Sham orbitals for each N, the 
arrows representing the sign of the spin. It can be seen 
that the spectra of both QDs exhibit the typical behavior 
found for conduction band electrons confined in axially 
symmetric QDs*i2ii£ with local maxima at shell filling 
values (N — 2,6) and at half- filling values (N = 4). 
These results are consistent with a sequential filling of the 
hole energy shell, according to the Aufbau principle and 
Hund's rule. It is worth noting that this result holds in 
spite of the fact that our model QDs are more correlated 
than the experimental ones. For instance, the Coulomb 
blockade energy between the two first holes in the ex- 
periment is 24 meV, similar to the value of 20.2 meV we 
obtain for the smaller QD. In contrast, the kinetic energy 
of our calculations is fairly smaller (about one third) than 
that inferred from the experimental data, in part owing 
to the assumption of pure heavy hole behavior. Even 
for the larger QD, no strongly spin-polarized system is 
predicted. At this point, one can notice that filling the 
hole energy shells sequentially means that the fifth and 
sixth holes occupy orbitals with the same absolute value 
of the envelope angular momentum z-projection as the 
third and fourth ones, namely \m z \ = 1. Therefore, the 
experimental observation that the magnetic field disper- 
sion of the iV=5 and JV=6-hole states is twice as strong 
as that of the 7V=3 and N=4 ones can no loger be ex- 
plained in terms of different angular momenta within the 
one-band model. However, different angular momenta 
may still be possible within a multi-band model, as HH- 
LH interaction yields hole states with different mixtures 
of envelope angular momenta, which has a dramatic ef- 
fect on the magnetic field dispersion^ We thus proceed 
to investigate role of HH-LH interaction in our model 
QDs. 

We use a three-dimensional, four-band k ■ p Hamilto- 
nian in the envelope function approximation, with a mag- 
netic field applied along the growth direction* 15 The QD 
we investigate has the same dimensions and confinement 
potential as the smaller dot of previous paragraphs. We 
use InGaAs Luttinger parameters (71=11. 01, 72=4.18, 
73 =4.84) which yield effective masses of 777^^=0.377 
and 777,2^=0.052. Figure 2 (solid lines) shows the lowest- 
lying single-particle energy levels as a function of the 
magnetic field. The levels are labeled by their total an- 
gular momentum z-projection, F z = m z + J z , where 
m z = 0, ±1, ±2... and J z = ±3/2, ±1/2 are the en- 
velope and Bloch angular momentum z-projections, re- 
spectively. The slopes of the energy levels vs. B are the 
result of the particular mixture of envelope angular mo- 
menta for each level (see Table I). For example, the pre- 
dominant component of the F z — 3/2 state at B = is 
that of J z = 3/2 (about 97.5% of the charge density), 
which is associated with m z = 0. This component is in- 
sensitive to the linear term of the magnetic field (orbital 



Zeeman splitting). 19 However, there is also a small con- 
tribution (2.5%) coming from the components associated 
with m z = 1, 2, 3, which introduces a slight dependence 
on the magnetic field. Likewise, for F z — 1/2 the lead- 
ing component (87%) is J z — 3/2, which is associated 
with 77i 2 = — 1. However, a non-neglegible J z = —1/2 
component (7%) with 777 z = 1 partially compensates for 
the m z = — 1 magnetic field dispersion, so that the final 
slope is not as steep as it would be in a one-band model. 
Finally, the F z = 5/2 state has a predominant (93%) 
J z = 3/2 component, connected with m z — 1, and minor 
components with m z = 2, 3, 4, which enhance the mag- 
netic field dispersion as compared to the one-band model. 
Similar reasonings can be made for F z = —3/2, —1/2 and 
—5/2 states. For a more accurate description of the mag- 
netic field dispersion, we need to add the spin Zeeman 
term, k \xb B J z , on the diagonal elements of our four- 
band Hamiltonian. In this term, \ib is the Bohr mag- 
neton and k the Lande factor for holes, which we take 
close to its value for GaAs, k = 1.2. 20 It is important to 
note that the spin Zeeman term, often neglegible in the 
description of conduction band electrons, may become 
comparable to the orbital Zeeman term for low-lying hole 
states. This is because such states have strong heavy hole 
character and thus a rather large effective mass acting 
upon the linear term of the magnetic field. The result 
of adding the spin Zeeman energy is illustrated in Fig- 
ure 2 (dashed lines). It can be seen that the magnetic 
field dispersion of the \F Z \ — 3/2 and \F Z \ = 5/2 states 
is slightly quenched, while that of the \F Z \ = 1/2 states 
is enhanced. The underlying reason is that the signs of 
J z and ?77 2 in the leading component of the \F Z \ = 1/2 
states are opposite. 

Next, we calculate the energy of the charging peaks 
vs. B which would follow from filling the energy levels of 
Figure 2 in an independent particle scheme, fJ,(N, B) = 
E(N, B) - E(N -1,B). To compare the different slopes 
of the individual charging peaks, we substract the peak 
energy at B — and represent the absolute value in Fig- 
ure 3. By comparing the figure to Figure 3 of Ref0 in the 
same range of magnetic field (t3=0-10 T), a fairly good 
agreement is observed. For example, the average slope of 
the first and second peaks is 0.04 meV/T (close to the ex- 
perimental value of 0.03 meV/T), the average slope of the 
third and fourth peaks is 0.18 meV/T (0.14 meV/T in the 
experiment), and the average slope of the fifth and sixth 
peaks is 0.32 meV/T (0.33 meV/T in the experiment). 
Our results therefore suggest that it is HH-LH interaction 
rather than Coulomb interaction the key factor on the ex- 
perimentally observed behavior, although of course the 
latter may help to improve the quantitative description 
of the system. Finally, we would like to point out that for 
QDs with other dimensions than those of the dot we have 
investigated here, the order of the excited single-particle 
energy levels may differ. However, agreement with the 
experimental results is only achieved for the energy level 
sequence shown in Figure 2. We would like to note 
though that a very recent theoretical work also investi- 
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gating the hole electronic configuration of self-assembled 
InAs QDs, albeit in absence of magnetic field, [L. He, G. 
Bester and A. Zunger, arXiv:cond-mat/0505330 points 
out that Coulomb interaction induces a change in the hole 
shell filling sequence. However, the QDs they study are 
significantly smaller than those expected from Refill and 
in any case the resulting electronic configuration upon in- 
clusion of the Coulomb term is equivalent to that we find 
in Figure 2 at an independent particle level. 
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FIG. 1: Addition energy vs. number of holes in a QD with 
R — 16 nm and h = 4 nm (solid line) and in a QD with 
R — 20 nm and h = 6 nm (dashed line) . The insets show the 
occupation of the Kohn-Sham orbitals for each N, the arrows 
representing the sign of the spin. 
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FIG. 2: Single-particle, four-band energy levels vs. magnetic 
field for a QD with R — 16 nm and h — 4 nm. Solid (dashed) 
lines are used for the energy levels calculated without (with) 
spin Zeeman splitting. 
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TABLE I: Envelope angular momentum z-projections and 
charge density weights associated with each J z component 
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tHe states with positive F z are shown. States with neg- 
have reversed order of weights and m z , and the sign 
c langes. 
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FIG. 3: Absolute value of the magnetic field-induced energy 
shift for charging peaks of 1 to 6 holes vs. the magnetic field. 
Solid lines represent the average slope of each pair of peaks. 



